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ON THE ELECTRON ENERGY DISTRIBUTION INDEX OF 
SWIFT GAMMA-RAY BURST AFTERGLOWS 
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ABSTRACT 

The electron energy distribution index, p, is a fundamental parameter of the synchrotron emission 
from a range of astronomical sources. Here we examine one such source of synchrotron emission, 
Gamma- Ray Burst afterglows observed by the Swift satellite. Within the framework of the blast wave 
model, we examine the constraints placed on the distribution of p by the observed X-ray spectral 
indices and parametrise the distribution. We find that the observed distribution of spectral indices 
are inconsistent with an underlying distribution of p composed of a single discrete value but consistent 
with a Gaussian distribution centred at p = 2.36 and having a width of 0.59. Furthermore, accepting 
that the underlying distribution is a Gaussian, we find the majority {> 94 percent) of GRB afterglows 
in our sample have cooling break frequencies less than the X-ray frequency. 

Subject headings: Gamma-ray burst: general — Radiation mechanisms: non-thermal — Acceleration 
of particles — Shock waves — Methods: statistical 
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1. INTRODUCTION 

The afterglow emission of Gamma- Ray Bursts (GRBs) 
is generally d escribed by the b last wave model 
()Rees fc Meszaro s 199^ iMeszaros et a l. 1998) which de- 
tails the temporal and spectral behaviour of the emis- 
sion that is created by external shocks when a colli- 
mated ultra-relativistic jet ploughs into the circumburst 
medium, driving a blast wave ahead of it. Fundamental 
to this model is the electron energy distribution index, 
p; a characteristic parameter of the process by which 
the electrons are accelerated to relativistic speeds and 
by which they radiate via synchrotron emission. This 
acceleration mechanism, common to many astronomical 
jet sources (as well as particle acceleration in the so- 
lar wind and supernovae, and the acceleration of cos- 
mic rays), is tho ught to be Fermi diffusive shock ac- 
celeration (iFermil Il 954) du e to the passage of an ex- 
ternal shock (iBlandford fc Ostriked Il978t iRieger et al.l 



|2007| ) after which the energy of the electrons, E, fol- 
lows a power-law distribution, dN (x E~PdE, with 
a cut-off at low energi es. This i s con sistent with 
recent PIC si mulations (iSpitko vskv 20 08|) and Monte 
Carlo models (lAchte rberg et al.! 2001: .Ellison fc Double! 
I2002t iLemoine fc Pe lleticr 2003) but at odds with others 
(jNiemiec fc Ostrowskiii2006c iNiemiec et al.l[2006[ ). 

The blast wave model describes how synchrotron emis- 
sion from relativistic electrons produces a smoothly- 
broken, broad band spectrum that is well characterised 
by a peak fiux and three, time evolving, break frequen- 
cies (peak frequency, Vn^; cooling frequency, i^c] syn- 
chrotron self-absorption frequency, v^ ) as well as the 
electron en ergy distribution index, p (jSari "eTalllTQM 
iGranot fc^S ari 2002). The spectrum is divided into 
four regimes by the three break frequencies and within 
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each regime the spectrum is asymptotically described by 
Fi, oc , where the spectral index, /?, is a function 
of p only. By comparing the observed X-ray spectra to 
the predicted asymptotic values of the synchrotron spec- 
tra, we can extract information about the electron en- 
ergy distribution index, p, which is dependent only on 
the underlying micro-physics of the acceleration process. 
Some (semi-) analytical calculations and simulations in- 
dicat e that there is a nearly universal value of ^ 2.2 — 2.4 
(e^ lKirk erall[2000l : lAchterberg etall[200ll : iSpitkovskvl 
120081) though other studies suggest that there is a la rge 
range of possible values for p of 1.5 — 4 ()Barins!jl2004[ ). 

Observationally, different methods have been applied 
to samples of BATSE, BeppoSAX and Swift bursts 
which reached the conclusion that the observed range 
of p values is not consistent with a single centra l 
value of p (Chevalier & Li' '2000*; Panaitescu & Kumai 
200 2: Shcn et al. 2006; Starhng ct al. 2008; Curran et a! 
20091) . The latter three showed that the width of 
the parent distribution is ~ 0.3 — 0.5. How- 
ever, in the studies so far there ha ve been some 
limitations: multi-wavelength st u dies (iChev alier fc"Lil 
2^00; Panaite scu fc Kumai) [20021: iStarling et al. 200^^ 
Curran et al.l 12009^) suffer from limited samples (< 10 
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sources each) with sufficient temporal and spectral ob- 
servations, while those studies that rely on X-ray after- 
glows alone are subject to a large uncertainty because the 
position of the cooling frequency, relative to the X-ray, 
is unknow n. The only X-ray study of Swift afterglows 
so far (Sh en et al.| [2006) used a very limited sample of 
spectral indices 30), dictated by the number of GRBs 
observed by Swift at the time, to estimate the distribu- 
tion of p. Neither, they did not take a statistical ap- 
proach to the position of the cooling frequency relative 
to the X-ray regime, as we do here. We interpret a much 
larger (^ 300) and, statistically, more significant sam- 
ple of Swift observed GRB afterglows, to constrain the 
distribution of the values of electron energy distribution 
index, p. In ^}2]we introduce our method while in fj3]we 
present the results of our Monte Carlo analyses and their 
implications in the overall context of GRB observations 
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and particle acceleration in general. We summarise our 
findings in 21 All uncertainties are quoted at the la 
confidence level. 

2. METHOD 

Our general method is to constrain the electron en- 
ergy distribution index, p, from the values of the X- 
ray s pectral indices, /3x or /3, observed by the Swift 
XRT (jBurrows et al.ll2005l) and detailed bv lEvans et"al] 
(|2009l Table 7) . A normalised histogram of the spec- 
tral indices of these 301 GRB spectra is plotted in fig- 
ure [TJ We derive p from the spectral index as op- 
posed to the temporal index because for a given spec- 
tral index, assuming the asymptotic limit, there are 
only two possible values of p depending on whether 
the cooling frequency, Vc is less than or greater than 
the X-ray regime, vx, while for a given temporal 
index there are multiple possible values which are 
model dependent (e.g. , the simple blast w ave model 
(|Rees fc Meszarosll99alM eszaros et al.lll998l ), or modifi- 
cations thereof (e.g.JGran ot & Kumar 2006; Genet et al.l 
[20071 IGhiselh ni et al.1 120071: lUhm fc Beloborodovfl2007l) 7 
for a discussio n on the choice of X-ray spectral index 
to derive p see ICurran et al.l 120091) . In accordance with 
synchrotron emission predicted by the blast wave model, 
we ascribe the behaviour of the unabsorbed X-ray spec- 
trum to be a single power law where the flux goes as: 
F^(v) oc and (i is the spectral index. Under the stan- 
dard assumptions of slow cooling and adiabatic expan- 
sion of the blast wave, the electron energy distribution 
index is related, in the asymptotic limit, to the spec- 
tral index by ei ther p = 20 [vc < vx) or p — 2(3 + 1 
ivc > vx) fe.g.. iGranot fc Sarill2002[ ): implying a differ- 
ence between the slopes of the two spectral regimes of 
A/3 = 0.5. Throughout we use the regime probability, 
X, as the probability that the cooling frequency is less 
than the frequency of the X-ray regime (i.e., Vc < ^'x) and 
1 — X is the probability that i/c > vx- We neglect the case 
where the cooling frequency may be passing through the 
X-ray regime (since there is no sign of spectral evolution 
in the sample) and the cases where the peak frequency, 
fc'in, or self absorption frequency, v^^ is greater than the 
X-ray regime as this is not observed in late time after- 
glows. 

To parameterise the underlying distribution of the elec- 
tron energy distribution index, p, from the X-ray spectral 
indices observed by Swift we use a maximum likelihood 
Monte Carlo method. This method uses a maximum 
likelihood fit to return the most likely parameters of the 
assumed underlying model, the errors on which are esti- 
mated via a Monte Carlo error analysis. Another Monte 
Carlo analysis tests the probability that the observed dis- 
tribution of spectral indices could be obtained from an 
underlying distribution of p described by the most likely 
parameters. 

In this method we first assume a model for the un- 
derlying distribution of p which we transform into a dis- 
tribution in spectral index space, via the regime proba- 
bility, X. We convolve this distribution with the mea- 
sured probability of the data set and calculate the log- 
likelihood of the parameters of the underlying model (see 
appendixE| . To estimate the most likely parameters of 
the underlying model and the regime probability, X, we 
minimise the log-likelihood (equates to maximising the 




Fig. 1. — Normalised histogram of the data (301 measurements) 
overlaid with high-resolution normalised histograms of the synthe- 
sized data sets (10^ X 301 data points) from the most likely param- 
eters of a single discrete p (grey line) and a Gaussian distribution 
of p (blue line) as detailed in Table[T| 

li kelihood) using the simulated annealing method (§ 10.9 
of I Press et al.|[l993 , and references therein). Uncertain- 
ties of the fit parameters are estimated via a Monte Carlo 
analysis, whereby the observed data are randomly per- 
turbed within their (asymmetric) Gaussian errors and 
refit multiple times (10^); the standard deviation of re- 
turned most likely parameters is used as a measure of 
the uncertainties. We then find the probability that the 
observed distribution of spectral indices could have orig- 
inated from an underlying distribution of p based on the 
most likely parameters estimated from the log-likelihood 
minimisation. We do so by generating 10** synthetic data 
sets drawn randomly from the underlying probability dis- 
tribution of p, randomly transformed into spectral index 
space via the regime probability, X, and further ran- 
domly perturbed within the (asymmetric) Gaussian er- 
rors of each of our original data points. Each of these 
synthetic data sets is fit as above and the value of the log- 
likelihoods recorded; one would expect that if the data 
are consistent with the underlying model, at the la level, 
the original log-likelihood value should fall within the la 
distribution of the synthetic values. With 10^ synthetic 
data sets we can measure the percentage to an accuracy 
of two decimal places and rule out chance agreements to 
the 4cr level; the ~ 2 x 10^ synthetic data sets required to 
rule out chance agreements at the 5a level was considered 
too costly, computationally. 

There are two underlying models, or hypotheses, re- 
garding the data that we want to test: that the observed 
distribution of spectral indices, /?, can be obtained from 
an underlying distribution of p composed of i) a sin- 
gle discrete value (SDp) and ii) a Gaussian distribution 
(GDp). Details of these models and their log-likelihoods 
are discussed in the Appendix, [Xj 

3. RESULTS AND DISCUSSION 

3.1. Distribution of p 

The results of our analysis are detailed in Table[T]which 
shows our most likely parameters for electron energy dis- 
tribution index, p, the related Gaussian scatter, cTp, the 
probability that the cooling frequency is less than the 
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TABLE 1 

Results of our likelihood fits of the observed spectral index distributions to a single discrete value (SDp) and a 

Gaussian distribution (GDp) of p. 

model p CT^ X I % 

~5Dp 2:25 = 093 ms 

(2.24±0.03) - (0.91±0.02) (1287±77) 

[2.253±0.011] - [0.93±0.02] [58±28] 100.00 

GDp 2.36 0.590 1.000 347 

(2.40±0.03) (0.600±0.007) (0.991±0.013) (478±25) 
[2.36ili0.05] [O.59O1I1O.OI2] [0.99ifc0.02] [382ifc24] 7.01 

Note. — The most likely values of electron energy distribution index, p, the related Gaussian scatter, o-p, the probability that the 
cooling frequency is less than the X-ray frequency (uc < i''k), X, and the log-likelihood of that fit, I. Values in brackets are the average 
and error values from the Monte Carlo error analysis of perturbed data sets, while those in square brackets are the average and standard 
deviations returned from the Monte Carlo analysis of synthesized data sets. % is the percentage of synthesized data sets with a better fit 
than the original. 



X-ray frequency {vc < i^x), X, and the log-hkelihood of 
that fit, Values in brackets are the average and error 
values from the Monte Carlo error analysis of perturbed 
data sets, while those in square brackets are the average 
and standard deviations returned from the Monte Carlo 
analysis of synthesized data sets. % is the percentage of 
synthesized data sets with a better fit than the original. 

Though our most likely discrete value of p ~ 2.25 
agrees well with the predicted univers al value of p ^ 
2.2 - 2.3 (e.g., ;K irk et al. 2000; Achtcrb eFT et al.ll2cioiD . 
it is comprehensively rejected by our tests; the hypoth- 
esis that the observed distribution of spectral indices, 
/3, can be obtained from an underlying distribution of 
p composed of a single discrete value (SDp) is rejected 
at the 4(7 level as all synthesized data sets had better 
log-likelihoods. The hypothesis that the observed dis- 
tribution can be obtained from an underlying Gaussian 
distribution (GDp) centred aX p = 2.36, having a width 
of 0.590 and regime probability, X — 1.000, is accept- 
able at the 1.5(7 level (Figure^]). As a visual aid, we 
compare the normalised histogram of the observed data 
with the high-resolution normalised, average histograms 
of the 10^ synthesized data sets (Figure[T|). Note that 
the SDp model is clearly a poor fit and exhibits a sec- 
ondary peak at /3 ~ 0.6 due to the fact that X = 0.96 for 
the likelihood fit of that model to the observed spectral 
indices. This result confirms the resu lts from previous 
small-sample GRB afterglow stu dies (S hen et al.l 120061 : 
iStarUng et all l2008t iCurran et al. 2009) as regards the 
non-universality of p, the central value at p ~ 2.0 — 2.5, 
and the width of the distribution of (7p ~ 0.3 — 0.5. How- 
ever, our results are based on a sample of bursts an or- 
der of magnitude larger than these studies and we took a 
statistical approach to the position of cooling frequency 
relative to the X-ray, by using the regime probability, X. 

Given that the value of p is not universal, it may be 
possible that it changes suddenly or evolves gradually 
with time or radius even in a single event as environ- 
mental shock parameters (e.g., magnetic field, ambient 
densi t y) change or evolve (e.g., iHamilton fc PetrosianI 
[l99l IVlahos et al.l[200l |Kaised|200l. It is also pos- 
sible that different components of a struct ured jet 
(jMeszaros et al l 119981: E umar fc Grano tI I2003I) mu lti- 
component jet (Berger et al. 2003; Huang et al.ll2004l ) or 
jet-cocoon (jZhang et al.l 12003 ) could have different val- 
ues of p. If either are the case our derived values of 
p should be considered time averaged values of the pa- 
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Fig. 2. — Histogram of the log-likelihood values of the 10* syn- 
thesized data sets generated from the most likely parameters, as- 
suming a Gaussian distribution of p, as detailed in Table[T] The 
log-likelihood of the original fit to the unperturbed data is marked 
along with the Icr distribution (dashed lines), showing that the 
data are consistent with the model at the 1.5a- level. 



rameter, though a change or evolution of p should be 
observable as a change or evolution of the synchrotron 
spectral index, /?, and no significant example of such an 
evolution has been observed in GRB afterglows. 

3.2. Limits on regime probability, X 

Though previous multi-ban d (opt ical - X-ray) stud- 
ies (e.g.,E naitcscujLKuinar 2002; Starling et al1l2008l: 
ICurran et al.. .2009) have shown that the cooling break 
frequencies of a number of GRB afterglows are greater 
than the X-ray frequency {u^ > i^), here, accepting that 
the underlying distribution of p is a Gaussian, we find 
that this number is consistent with zero {X = 1). From 
our Monte Carlo error analysis (where we refit the ran- 
domly perturbed data multiple times) we can plot the 
distribution of possible values of X (Figure|3J note that 
this is a log-log plot). The distribution of the errors 
is clearly not Gaussian so we should not use the stan- 
dard deviation as a measure of error as we have in Ta- 
ble[Tl We can however place a nominal 3cr lower limit 
on the value of X by estimating the value at which 
there is 99.7 percent containment; we find that this is 
at 0.935 < X < 0.936, compared to the nominal Icr 
(68.2 percent containment) at 0.988 < X < 0.989. Any 



4 



Curran et al. 




Fig. 3. — Histogram (on a log-log plot) of the values of X re- 
turned by the Monte Carlo error analysis, assuming a Gaussian 
distribution of p. The dashed vertical line shows the 99.7% con- 
tainment at X = 0.935 which we take as a measure of the lower 
limit on X. 



value of X less than this limit is inconsistent with the 
distribution of the data and would produce a clear sec- 
ondary peak at P — {p~ l)/2, not observed in the data 
(Figure[T|). To avoid this peak, an extremely wide dis- 
tribution is needed, much wider than the spread of the 
observed data. 

Hence, the upper limit on the percentage of GRBs in 
our sample where the cooling frequency is greater than 
the X-ray frequency {vc > i^x) is approximately 6.5 per- 
cent, or ~ 20 out of 301 GRBs. This is a discrepancy 
from the ratio observe d in t he multi-band studies (e.g., 
iStarling et al.l : fCurran et al.l : 5 out of 10 and 2 out of 6, 
respectively) but not seriously so given the extremely low 
number statistics of those studies. If we create a sam- 
ple of 6 random bursts from our statistical distribution of 
301, there would be a non-negligible (~ 5 percent) chance 



that 2, or more, bursts would have a cooling frequency 
greater than the X-ray freq uency, consist ent with the 
afore mentioned S wift study (jCurran et al.| ). The study 
of Sta rling et al.l is based on a sample of BeppoSAX, as 
opposed to Swift, GRBs which may have a different limit 
on the regime probability, X . An investigation as to why 
the cooling frequency of afterglows follow this apparent 
limit, or its implications regarding the distribution of 
other parameters, is beyond the scope of this work. 

4. CONCLUSION 

We use the X-ray spectral indices of gamma-ray burst 
afterglows observed by the XRT aboard Swift to param- 
eterise the underlying distribution of the electron energy 
distribution index, p, within the framework of the blast 
wave model. The electron energy distribution index is a 
fundamental parameter of the synchrotron emission from 
a range of astronomical sources and in this case of the 
synchrotron emission of GRB afterglows. We use a max- 
imum likelihood Monte Carlo analysis to test two hy- 
potheses, namely that the observed distribution of spec- 
tral indices, /?, can be obtained from an underlying dis- 
tribution of p composed of i) a single discrete value and 
ii) a Gaussian distribution. We find that the observed 
distribution of spectral indices are inconsistent with the 
first hypothesis but consistent with the second, a Gaus- 
sian distribution centred at p = 2.36 and having a width 
of 0.59. Furthermore, if we accept that the underlying 
distribution is a Gaussian, the majority {> 94 percent) 
of GRB afterglows in our sample have a cooling break 
frequency less than the X-ray frequency. 
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APPENDIX 
PROBABILITY FUNCTIONS 
The probability of spectral index, /3, given a measured value, Pi, with asymmetric errors, ai± (IStarling et al.l[2008l) 
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P^(/3|A,^»±) = 
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where ai± = [a.i^ + ai+)l2 and ai± 



Assume that the distribution of the electron energy distribution index, p, can be described by a Gaussian probability: 



P(pb,crp) = 
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27r (T, 



■ exp 
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2ol 



(A2) 



where p is the central, or mean, value of the distribution and dp is the standard deviation of the distribution. The 
distribution of the observed X-ray spectral indices, /?, can then be described by a double Gaussian: 



P(/3|p,ap,X) = 



1 



2tx ((Tp/2) 



X exp 



-(/3-p/2)2 



2(ap/2)2 



(1 - X) exp 



-(/3~p/2 + 0.5)^ 
2(ap/2)2 



(A3) 

where X is the probability that the the cooling frquency is less than the X-ray frequency (vc < z^x), and (\ — X) is 
the probability that the cooling frequency is greater than the X-ray frequency (y^ > i^x)- Convolving this probability 
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with the measured probability (equation IA1[) gives: 

1 1 



/2tt 
X I X exp 



y^(ap/2)2 + ai±2 

-(A-p/2)2 



.2((V2)2 + af±) 
(Ji- {(3i < p/2; first exponential) 



+ {l-X) exp 



-(A-p/2 + 0.5)^ 
L 2((V2)2 + a2^) 



(A4) 



where cr + = > P/2; first exponential) 

" ' (Ji- < p/2 — 0.5; second exponential) 
0"!+ > p/2 — 0.5; second exponential) . 
This leads to 3 different permutations of equation IA41 since the combination of Pi < p/2 — 0.5 and /3i > p/2 does not 
exist. The log-likelihood of these probability distributions is given by 

l{p,ap,X) = -2\nL{p,ap,X) 

N 

^~2\nl[P,{l3\p,ap,X,l3„a,±) 

i=l 

N 

= -2j2 In-P»(/3|P, 'Jp, X, 13,, a,±) 
1=1 

N 

= N ln(2^) + J2 ln(('^p/2)' + ^') 

i=l 

N , 

-2^ In ( Xexp 



-(/3.-p/2)^ 



2(K/2)2 + af±) 



(f - X) exp 



-(A -p/2 + 0.5) 



21 



2(K/2)2 + a2^) 



(A5) 

where L(p, ap, X) is the likelihood function. This is a numerically calculable function that is minimised to find the 
most likely parameters, errors on which can be estimated via a Monte Carlo error analysis. If the distribution of the 
electron energy distribution index, p, can be described by a single discrete value, 



P(p|p) 



_ / 1 ip = p) 
(p^p) 



, then the probability of /3 is 



(X {P^p/2) 
V{m^{ 1-^ (/? = p/2-0.5) 

i (/? ^ p/2 & /3 ^ p/2 - 0.5) 

and the above convolved probabilities and likelihoods hold with Cp = 0. 



(A6) 
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